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MADCAP IV 


A Multiplexed ADC and Analog Plotter Program 

Author: G. W. Dulaney- 

Digital Equipment Corp., Maynard, Mass. 

ABSTRACT: 

The program was designed as a foundation from which 
the construction of routines for data acquisition and 
analysis can be more easily done. By the use of short 
overlays, it is able to accept many forms of input data 
including analog output from peripheral equipment such 
as a mass spectrometer, ESR, NMR, curve follower, etc. 

This output is converted into digital form using the AX08 
Lab Peripheral's 9-bit ADC. The mode of input for mass 
spec, ESR, and NMR requires signal averaging to remove 
noise. This is accomplished by performing 'boxcar inte¬ 
gration,' which consists of averaging a certain number of 
consecutive samples to generate a data point in a controlled 
period of time. The lowest S/N increase for 60-cycle noise 
and one scan is about a factor of ten (128 samples per 
point at 50 points per second), and it can be as high as 
40 (2048 samples per point) for a slow scan of more than 
2.0 seconds per point. 
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Important, also, is the program's utilization of a 
general-purpose X-Y analog plotter in parallel with the 
AX08 oscilloscope output as opposed to the more expensive 
incremental plotter normally employed for computer output 
plotting. 

Output is in the form of printed (and/or punched) in¬ 
formation and a scaled plot which can be either a line plot 
or point plot with optional axes and tic marks drawn in. 

The MADCAP IV Basic Package accepts digital data input 
from paper tape or the keyboard for plotting. Included 
in the Basic Package are routines to compute the best 
line between adjacent points for plotting to prevent plot¬ 
ter servo overshoot and a one-page program modifier routine 
which operates much like ODT-8 in that locations in core 
can be changed to suit a special purpose. Without the modi¬ 
fier, about 180010 locations are available for overlays 
and data. 

Analog input routines are supplied as separate, short 
binary tapes to be added to the Basic Package. Overlays 
currently available include a flexible signal averager, 
a routine to use a photoelectric curve follower, and a 
specialized routine to accept data from a mass spectrometer 
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to generate ionization efficiency curves. Integration 
and curve smoothing routines also are available to massage 
data after acquisition. All data acquisition and reduc¬ 
tion routines are re-executable without reloading the 
binary tape. 

The Basic Package has been used as the basis for the 
LOPGAP Package which can input data, compute Lorentzian 
peak spectra, compute residuals between original and cal¬ 
culated data, integrate either type of data, strip out un¬ 
desired peaks, remove nonlinear baseline, and smooth data 
using an eleven-point parabolic function. Results of any 
of these operations can be displayed on the oscilloscope 
and/or plotted and are easily executed from the ASR key¬ 
board. LORCAP can be converted by means of a short overlay 
to handle Gaussian peaks. 

The LORCAP Package has been used in the analysis of 
Mossbauer and NMR data and with the addition of a Gaussian 
line shape overlay, it can find application in the simula¬ 
tion of IR, visible, and other forms of spectroscopic data. 

The author encourages response from users of the 
program with respect to how well it meets their applica¬ 
tions, any new routines they may wish to make available 
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for use with it, their satisfaction (or lack of it, too) 
with the concepts used, and any changes they wish to pro¬ 
pose. 

I. REQUIREMENTS: 

The system on which the program was written and used 
was a basic 4K PDP-8/I LAB-8 Computerpack with high-speed 
reader. Peripherals necessary (See Appendix 2) are the 
AX08 Lab Peripheral, ASR—33, and an oscilloscope. A 
Hewlett-Packard 2D-2 X-Y analog plotter with remote relay 
actuated pen lift and Y-axis retransmitting potentiometer 
and a Hewlett-Packard F3B photoelectric curve follower also 
were used with good results. A felt-tip pen was found 
advantageous in making point plots and, in fact, a simple 
sheet metal holder was fashioned which made use quite 
satisfactorily of dime-store refill cartridges. A relay 
driver and relay were added to the AX08 and triggered by 
S0 for pen up-down control on the plotter. 

Hardware desired and/or necessary for signal averaging 
are preamps with a gain of approximately 20 for accepting 
approximately 50 mv signals - a FLIP-CHIP type A200 opera¬ 
tional amplifier is suitable - and a relay output - FLIP-CHIP 
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type W800 is useful here. The scans are initiated by an 
external trigger (e.g., +1 to -3 volts) impressed on SI 
of the AX08. This can be done a variety of ways and is, 
of course, dependent on the nature of the signal source 
device. For NMR or ESR instruments unequipped for computer 
output, a microswitch actuated by the recorder carriage 
can be used to switch in a simple battery circuit to 
trigger the computer. The end-of-scan marker from S0 of 

/C 

AX08 (-3 to 0 volts for 1 sec., and then back to -3 V) 
can be divided down and inserted into the Y axis of the 
instrument's recorder to allow accurate calibration. 

These are merely suggestions for thought, and the user 
can obviously apply his ingenuity to better implement and 
utilize the signals. 

The arithmetic and I/O routines utilize Floating Point 
Package #2 (DIGITAL 8-5B-S), a relocated version of DECUS 
No. 8-44 Modified Output Controller, and overlays written 
by the author to allow use of a high-speed reader punch 
under control of C(56). The resulting FPP is herein 
called the Modified Fl. Pt. Pkg. #2 and has been incorpor¬ 
ated into both the Basic Package and the LORCAP Package. 
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II. USAGE: 

The FPP accepts all standard input formats; e.g., 
any of the following are accepted as 1.0: "1", "1.0", 

"1E0", "1.0E0", "0.1E1", and "10E-1". 

In each case, the number is terminated by anything 
other than a digit, "E", or a period. A space is suggested 
as a terminator. Leading spaces and spaces following a 
plus or minus sign are ignored, a RETURN not immediately 
preceded by a digit is ignored, and plus 

signs are converted to spaces on output. On input, a 
RETURN is not echoed with a LINEFEED. 

A number can be cancelled on input if RUBOUT is 
struck before a terminator is entered. An exclamation 
point is echoed;then the user can enter the correct value. 

The subroutine which handles the response to questions 
will accept only "Y" or "N" as valid input. All other 
characters generate "?" and are ignored. No terminator 
is needed. 

Insure that the BIN loader is in core - see CONSOLE 
MANUAL for instructions to load it. For all systems ex¬ 
cept LORCAP, load the Basic Package (which includes the 
Modified FPP) into desired reader, turn Teletype on line. 
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LOAD ADDRESS 7,77, set bit 0 = 0 to use high-speed reader 
and press START. The tape has only one checksum. Now 
load any desired overlays; e.g., the mass spec routine 
and/or the integrator, etc., by placing them in desired 
reader and hitting CONTINUE. 

The LORCAP Package xs a two-part binary tape, the 
second of which is a set of operating instructions that 
will be typed out when the program is started. The first 
section contains all that is needed to operate and the 
second section is optional. 

All data input routines, including LORCAP, use SA=0200. 
A heading is typed then the user can type in a title of 
any desired length, terminating it with ALTMODE. The dis¬ 
cussion is now separated into parts descriptive of each 
mode of operation. 

A. Digital data input: 

Answer "N" to ANALOG INPUT? and give the number of 
points to be retained for NO. POINTS=. Up to 500 points 
with X data or up to 1000 without X data can be entered. 

If X data is to be entered, answer "Y" to HAVE X DATA? The 
program then requests PTS. TO SKIP=, i.e., the number of 


points to ignore on read in. 


This feature allows the user 
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to select any contiguous section or block of points from 
the data tape. Scaling variables for Y data (and X, if 
present) are entered next. RANGE should be chosen to be 
at least as large as the difference between the largest 
and smallest values within the data block. MIN is the 
smallest value, or less, in the block. These are used to 
scale the data to fit within a 0-500 range for display 
and can be chosen to maintain exact intensity (and X 
axis) calibration, if necessary. If no X data is to be 
entered, X RANGE and MIN are not asked for. Then the 
user is asked TTY I/O? Load the data tape and prepare 
the desired reader before answering. "N" will start 
reading from the high-speed reader. When input is complete, 
scope display commences. See Section II.G for further 
instructions. 

B. Mass Spec Analog Input: 

This mode was designed for measuring relative ion cur¬ 
rent as a function of accelerating voltage, wherein the 
accelerating voltage is manually incremented on a pot. A 
voltage proportional to the ion current is sampled and 
averaged by injecting an external trigger pulse (+1.5 
going to -1.5 volts) into SI of the AX08 to begin each 
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data point. 2048 samples are averaged per point, and the 
results for each point can be printed out as soon as they 
are complete. This requires about 1.5 second per point. 

The currently accumulated data is displayed on the scope 
while waiting for the next trigger pulse. 

To enter this mode, answer "Y" to ANALOG INPUT and 
enter the desired MULTIPLEX CHANNEL. To print and/or 
punch the data as it is taken, answer "Y" to OUTPUT? and 
prepare the punch before answering TTY I/O?; "N" assumes 
the high-speed punch is to be used. In this mode, the 
number of points is left unspecified, and the user can 
enter up to 1000 points or terminate entry when desired. 

The program will punch leader tape if output was desired 
and then await the first trigger pulse. If output was not 
requested, the Teletype bell is rung when each point is 
complete. 

When it is desired to terminate data entry, set BIT 
11=1. Trailer is punched if output was requested, then the 
program types out the number of points taken before scope 
display commences. See Section II.G. 

C. Curve Follower Input: 

Note: For plotters other than HP 2D-2 and/or curve 

followers other than HP F3B, refer to manuals supplied by 


the manufacturer. 
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Connect the Y output of the F3B to Y input of recorder, 
place recorder in LF mode and connect the multipin plug 
into receptacle on rear of recorder marked RETRANSMITTING 
SLIDEWIRE. Make signal connection to AX08 and insert 
pickup head in pen holder. See manuals for more detail. 

Answer "Y" to ANALOG INPUT?, give desired MULTIPLEX 
CHANNEL, and select desired PT. INDEX (Table 1). An 
illegal PT. INDEX results in it being requested again. 

The program now types an instruction to set up X ZERO and 
GAIN on recorder to encompass desired region to be sampled. 
See Section II.G under PLOT for detail on the plotter cal¬ 
ibration routine. Axes can be drawn on the original, if 
desired. 

When range is set, turn RC clock COARSE to first or 
second position from slowest ( 0.5 to 1 sec/cycle) and type 
"0". Center pickup head over line, increasing follower 
SENSITIVITY and GAIN just short of oscillation. Type "G" 
to begin data input. When input is complete, the scope 
display routines are entered. 

Turn recorder to standby and set follower gain sensi¬ 
tivity fully counter-clockwise. Refer to Section II.G for 


output details. 
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D. Slow Scan Signal Averaging Input: 

In this mode, a remote trigger (+ TO -) impressed on 
SI of AX08 commences the scan after a specified delay time. 
Up to about 500 scans can be taken and averaged. After 
each scan, an outgoing trigger is available from S0 of the 
AX08. It is a -3 to 0 volt pulse lasting about 1 second. 

A new scan must again be triggered thru Si. 

Digital data input via Section II.A is unavailable 
with the Averager, and Y axis tic marks are not made. 

After starting at 0200, type a title, hit ALTMODE, 
and supply desired multiplexor channel, point index, and 
SEC/SCAN. 

If the input for PT. INDEX and SEC/SCAN give a value 
less than.016 sec/pt, new input is requested, starting 
from PT. INDEX. Then set the MAX SCANS «500), and, 
finally, the delay time in secs. (<400) (See Tables 1 and 2.) 

Input commences with the first trigger and continues 
to the end of the scan. While in progress, the scope will 
follow the ADC 1 s as they are performed with the X axis 
scaled to the number of points. 
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The run can be interrupted at the end of the scan 
in progress by setting BIT 0 of the SR to 1 at which time 
the number of scans taken is printed and scope display 
of stored data commences. The size of the display can be 
changed by typing "D" on the ASR to decrease size or "M" 
for magnification. Typing "S" to stop display allows 
the user the option of returning to accept more data (up 
to MAX. SCANS) or terminate the run as it stands by giving 
the desired answer to MORE SCANS? After termination, the 
scope routines are entered. If the integrator and/or one 
of the smoothers is desired, this can now be loaded with 
the BIN loader. Restart program at 0201 to reenter scope 
routines. See Section II.G for further details. 
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Computer memory cycle time is used for the time base; 
thus, accuracy will vary from one machine to another but 
reproducibility is quite good. 

Accuracy can be improved to about 0.5% using the 
following calibration method: 

1. Accurately measure the scan time, T, with 500 

points (PT.INDEX=3), 180 seconds or more for 

SEC/SCAN (defines ST below), and 0 for DELAY as 
the input. Since there is a one-second delay 
loop used for the end of scan marker, carefully 
watch the PC and MA lights to detect the actual 
end of the scan. For convenience, C(3605) can 
be changed from 6331(XRIN) to 7604(LAS) to allow 
SR bit 1 to trigger the scan. 

2. Compute new memory cycle time, MCT: 

MCT = (1.450X10 -6 ) xST/ t 
( 1.450 microseconds was the value assembled into 
the Averager.) 

3. A short routine is available at startup only by 
SA 0400 to allow resetting the MCT to a new value. 
The routine is destroyed by data. After entry 

of the calibration value, the routine exits to 


0200; i.e., the normal starting address. 
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PT. 


Table 1: 



Number of 

Points Assigned 



According 

to PT. INDEX 


INDEX 

Follower 

Averager Min 

.Sec./Scan 

1 

125 

50 

0.8 

2 

250 

100 

1.6 

3 

500 

500 

8.0 

4 


1000 

16.0 


Table 2: 

Number of ADC Samples per Point as 
a Function of Time Density for Averager; 


Seconds/Point 

Samples/Point 

.015 >0.5 

128 

0.5 > 1.0 

512 

1.0 > 2.0 

1024 


2.0 > 


2047 
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E. MODIFIER Routine (Resides in the Basic Package): 

The routine is similar to "ODT-8" but is used only 
for making changes in core locations and uses only one 
page of core. LA 2400, hit START, and the routine will 
give a carriage return-linefeed. To open a location, 
type the octal value then a space. The contents are then 
typed out followed by a colon. Insert the correction then 
type "C" to enter the correction. If no correction is 
desired, hit the SPACE bar or RETURN. To inspect the 
next location after closing one, type "N" and the location 
and contents are printed with a colon. To reopen the same; 
i.e., last, location used, type "S".RUBOUT will terminate 
any operation and restart the routine. 

A jump into any part of core can be made by typing 
"R" for return, then the location and a SPACE. (Be sure 
to set any controlling SR bits used in that routine.) 

Of the currently available routines for use with 
MADCAP IV, only a 1000-point signal averaging run and 
LORCAP destroy the MODIFIER. It can be reloaded after 
averager data acquisition, but is never available to LORCAP 
in its normal configuration. 
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As may often happen in certain experiments, a data 
set may contain a few obviously "bad" points. If allowed 
to remain, they will mar performance of the scaling, 
smoothing, and other routines. Roughly tagging their po¬ 
sitions on the scope, they can be located and reset using 
the MODIFIER. Unless changed, the data array starts at 
location 400g. 

F. Utility Routines for MADCAP IV: 

1. SMOT29 and SMOTll: 

Two data smoothing routines are available which can 
be loaded at start time with all but a 1000-point signal 
averaging run (load after the run) or LORCAP (already 
contains SMOTll). Both are modifications of M. W. King's 
LESQ29 and LESQ11 (DECUS No. 5/8-69). Only one or the 
other may be present at a time. 

To execute smoothing type FILTER: (or F:) during 
scope display. Up to 50 seconds may be required for 1000 
points with SMOT29. When the operation is complete, the 
scope routines are entered. The data can be smoothed as 


many times as desired. 
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2. INTEG: 

This is a numerical integration routine using the 
second order Simpson's Rule. It is executed by typing 
INTEGRATE: (or I:) during scope display. It requests a 
value for a PIVOT LINE with a value in a range of 0-500. 
This value will be subtracted from each data point. In 
essence, the routine computes the running area around the 
pivot at each point then scales the data down to a range 
of 500. The factor used is printed out then the scope 
routines are entered. It can be executed as often as 
desired, each time supplying the appropriate PIVOT LINE. 
Like the smoothers, a 1000-point signal averaging run 
will destroy the integrator; but it, too, can be loaded 
after the run. It requires about 5 seconds for 1000 
points. LORCAP contains its own version of INTEG. 

G. Program Control During Scope Display: 

All routines of MADCAP exit to the scope display 
routines. During the display, the Teletype keyboard is 
active and will accept a variety of commands which initiate 
data treatment, output, and new input. Hence, once the 
Basic Package is started, it becomes independent of the 


switch register and is totally keyboard oriented. Certain 
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of the analog data input overlays do, however, require 
use of the switch register. 

The Basic Package allows the following commands. 

Each begins with a unique letter but allows any number of 
characters to be entered after that character until receipt 
of a colon which executes the command, if it is legal, or 
a RUBOUT which deletes the command. Illegal; i.e., unde¬ 
fined, commands and RUBOUT force a return to the scope 


display. 


SCALE; 

Scales the Y data to a range of 0-500(10) 
and prints out the scaling parameters 
used. 

OUTPUT: 

Outputs the Y data onto the device to 
be specified by the user in answer to 
TTY I/O?. "N" assumes use of the 

high-speed punch. In either case, 
prepare the device (turn it on, run 
leader, etc.) before answering the 
question since output commences im¬ 
mediately upon receipt of "Y" or "N". 

PLOT: 

Plots the data using a solid line or 
as points depending upon the response 
to LINE PLOT?. "N" assumes a point 
plot is desired. Upon receipt of "Y" 
or "N", the plotter calibration 
routines are entered. In this sub¬ 
mode, the keyboard is again used to 
direct control. The legal characters 
are X, Y, 0, M, and G and are typed 
without colons, etc., simply as the 
single character. Before typing M or 

G, set the RC clock COARSE switch on 
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PLOT: 
(cont'd.) 


first or second from slowest position. 
This controls pen speed. During draw¬ 
ing of lines, the COARSE switch can be 
changed, but do not move it during 
point plots or when drawing tic marks. 
The FINE control potentiometer can 
be changed anytime. Because of the 
0 to -10 volt display signals, it is 
convenient to use the following cali¬ 
bration sequence. Type "X" (gives 
approximately 0 volts) and set the 
X zero pot of the recorder; type "0" 
and set the X variable range pot; 
repeat until satisfactory then cali¬ 
brate the Y axis similarly. 

X = Move pen to maximum X coordi¬ 
nate . 

Y = Move pen to maximum Y coordi¬ 
nate. 

0 = Move pen to X-Y origin. 

M = Mark off a frame to inscribe 
the plot and put quadrant 
tic marks on the X and Y axes. 

G = Exit the calibration routine 
and commence plotting. 

When the plot is complete, a message 
is printed requesting that the plotter 
be turned off (or put on standby). 

Do so immediately, since the scope 
display is restarted which will drive 
the pen wild in its effort to follow 
the scope analog signals. The key¬ 
board is live during plotting, so un¬ 
less premature termination is desired 
(shut plotter off then hit SPACE bar 
and a colon to do so), do not depress 
any of the keys. A simple way to arrest 
the scope display after the plot is to 
hit the SPACE bar as soon as it starts, 
shut plotter off, then type a colon. 
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REENTER: 


ANALOG: 


Allows reentry of digital data; i.e., 
it jumps to the routine to allow the 
user to enter a new set of data points 
from tape or Teletype. Note: This 
option is not available to the Averager. 

Restarts an analog data input overlay, 
if present, and if not, the command is 
equivalent to REENTER. 


The integrating and smoothing overlays, if present, 
can also be initiated. See Section F. If they are not 
present, the scope display is restarted. 


FILTER: Filters or smooths the Y data using 

SM0T11 or SMOT29 (11- or 29-point para¬ 
bolic least squares fitting). 

INTEGRATE: Integrates the Y data after subtract¬ 

ing the user supplied value for PIVOT 
LINE, scales the integral to a range 
of 500(10), and prints the scale down 
multiplier. 

In all cases, only the first character of a command 
is decoded; thus, the single letter followed by a colon 
will suffice. Exercise caution in doing so. 

If/ for some reason, the user must stop the program, 
depress the STOP console switch. The scope display can 
be restarted by SA=0201 or a program restart can be made 
by SA=0200. 
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H * LORCAP IV Package: 

1. General Description: 

As described in the abstract, LORCAP is a set of 
data processing routines as opposed to data acceptance. 

It was written to work with Lorentzian spectra as can be 
obtained from Mossbauer and NMR spectrometry and can 
be a useful model to simulate IR and UV data. The func¬ 
tional Lorentzian expression relating height, width, and 

arbitrary X-axis position is as follows: 

N 

Yi = HI HjWj 2 /( (Xi— Pj)2 + Wj 2 ) 

j 

where Yi is the sum of intensities at point i over N peaks 
characterized by their independent Hj 's, Wj 1 s, and Pj's. 

The result for a single peak is an envelope similar to a 
Gaussian peak but having more extended "tails". The 
width, defined as the half—width at half—height, corresponds 
to a Gaussian sigma value, the height and position have 
identical correspondence, and the area is directly propor¬ 
tional to the height-width product. 

LORCAP will handle up to 1000 points total including 
X data, if needed. Thus, up to 500 but an equal number of 
points for each type (raw and calculated) or 1000 points 
of either type, all without X data, may be entered. 
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If X data is desired, it requires on array equal to the 
number of points and the spectrum calculator will use it 
if its presence was so declared in answer to HAVE X DATA?. 

If the answer was negative, an equally incremented X axis 
scaled to X RNG and XI is used. 

2 . Usage of LORCAP: 

The binary tape is in two parts. Load the second 
only if the operating instructions are desired. This 
section was included to aid those users who, like the 
author, usually forget, loan, or otherwise incapacitate 
the manual or book needed while on the system. The instruc¬ 
tions are loaded into the data array and are destroyed 
when data is entered. Hence, they are typed out only when 
the program is first started. Similarly, but within the 
first section, a short dialog subroutine, which loads 
into the peak parameter array and allows several structural 
changes to be made, is used only at startup. 

LORCAP uses SA 0200, types a heading, then asks for 
the initial data block number. Two data blocks are assumed, 
0 for raw data and 1 for calculated or LORCAP data. If raw 
data is to be entered first, type 0, or type 1 if a spectrum 
is to be calculated first. 
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The program then desires some information which is 
used to change the structure of LORCAP slightly. The 
first choice allows the spectrum calculator to compute 
absorption dips (descending from a baseline) or peaks 
(growing up from baseline). Supply the desired answer 
to WANT ABSORPTION 'DIPS'?. The next option allows the 
assignment of equal widths to all peaks, hence, requiring 
the width be entered only for the first peak when the 
parameters are entered. Answer ALL WIDTHS EQUAL? as de¬ 
sired. If "N", a width must be entered for each peak. 

Both options are fixed until the program is reloaded. If 
the second part of the binary tape was loaded, a brief 
set of LORCAP's operating instructions will now be printed 
out. 

After the initial dialogue, CALC'N INPUT? is asked. 

To enter raw data, answer "N". Raw data entry is identical 
to that described in Section II.A. If "Y", a set of 
variables and parameters are to be entered and used to 
calculate a spectrum. Give the NO. POINTS to be assigned 
to each array, answer HAVE X DATA? to describe the raw 
data or give "N" if none is to be entered, and type in the 
NO. PEAKS to be convoluted (up to 25). The X RNG, XI, INT.MPLR, 
and BASE are then requested. The X axis used in the calculation 
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can be scaled to any desired units by appropriately set¬ 
ting the range (difference of right from left extreme) and 
first X value (left extreme) for X RNG and XI. If X data 
is present, use the same values given for X RANGE and XMIN. 
The INT.MPLR is an arbitrary multiplier of all heights 
and the baseline. This allows scaling the Y axis to 
suit. The product of INT.MPLR and the sum of the largest 
height and BASE must be less than 500 to prevent scope 
display overlap. For example, to allow entry of heights 
scaled to a largest peak of 100 with a base of 10, enter 
a value of 4.555 (500/110) for INT.MPLR to get a fully 
scaled display, assuming no intensity is acquired by over¬ 
lap of peaks. 

Finally, the table of heights, widths, and positions 
for all peaks is entered. If the user specified all 
widths to be equal, width is entered for only the first 
peak, and then HGT followed by POS for all subsequent 
peaks. When the parameters for a peak are complete, a 
carriage return—linefeed is issued, and when all input is 
complete, the message COOL IT. is typed. The program is 
now computing the spectrum and will require about 0.06 
seconds per point per peak; e.g., five peaks and 1000 
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points take about five minutes to compute. When complete, 
a carriage return-line feed is issued and scope display 
commences. 

During scope display, LORCAP uses an expanded version 
of the Basic Package keyboard command decoder. Each 
command begins with a unique pair of characters as opposed 
to one and block 0 and 1 can be alternately displayed by 
typing "0" and "1". Illegal commands generate "?" and 
are ignored. As in the Basic Package, the command is not 
tested for validity until a colon is received and a RUBOUT 
will cancel it and return to the scope display. 

The following four commands assume that raw data is 
in block 0 and calculated data is or will be in block 1. 
Therefore, the user should be displaying block 1 before 
entering one of these four commands (even though no real 
data may be in that array). 

COMPUTE: Computes a spectrum from a new set of var¬ 

iables and parameters to be entered. This 
is logically equivalent to answering "Y" 
to CALC'N INPUT? on startup. 

ALTER: Alters the parameters indexed by their 

PEAK & PARA #'S where the peak number is 
indexed in the order of initial entry 
(it is not a function of the parameters), 
and the para numbers are as follows, even 
if the equal widths modification is in 
effect; i.e., although entered as equal. 
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the widths of all peaks are handled inde¬ 
pendently and can be reset to any desired 
value. 

N, 1: height of peak N; 

N, 2: width of peak N; 

N, 3: position of peak N. 

Also 

0> 1: intensity multiplier; 

O, 2; baseline. 

Current values for the parameters can be 
printed out thusly: 

N, 0: prints parameters for peak N; 

O, 0: prints out all parameters. Exit 

is made without computing the 
spectrum in either case. 

When desired alterations are complete, type: 

0/~lt exit alteration mode and compute 
a spectrum from the current para¬ 
meters . 

SWAP: Swaps block 0 and block 1 data. 

SUBTRACT: Gives point-by-point subtraction of block 

1 data from block 0 and adds 250 1O to each 
resultant point. Double subtraction will 
restore the original array. 


^■11 the following operations can be performed on 
whichever array is currently displayed. For brevity, CDA 
is defined as the currently displayed array. In all cases, 
including the above four commands, the program will use the 
latest value supplied for the block number; i.e., if it 
has been already set to the desired value previously, it need 
not be reinserted, except in the interest of safety. 




STRIP: 


On entering the first time, a pair of 
bright spots are placed on the scope. 

Their positions are controlled by analog 
input potentiometers through multiplex 
channels 34-37 of the AX08. The top pair 
control the "left side" point's X and Y 
coordinates, resp., and the bottom pair 
similarly govern the "right side" point. 
Three options will be available: (1) 

replace the data with the best straight 
line from one point to the other; e.g., 
strip out a peak; (2) subtract a sloping 
baseline from all data points; and 
(3) subtract the straight line between 
the points from the data between points; 
i.e., partial baseline restoration. 

For option one, carefully place the two 
test points to remove the desired section 
of the data. They will normally be some¬ 
where on the data curve. The X positions 
need not be precise for option two. 

Simply place the test points with the 
proper slope and amplitude, usually lying 
just under the data curve. The points 
do not have to be located at the X axis 
extrema. Data processing using option 
three often involves performing option 
two first and then removing the residual 
non-linear baseline in short segments via 
option three. Position the points as 
would be done for option one, but usually, 
the points should not lie on the data but 
be somewhat below it. One, in fact, may 
often lie near the bottom of the screen 
as in the case of NMR data treatment, for 
example. After the points are positioned, 
the STRIP routines must be reentered, 
this time to select and execute one of 
the above options. To execute option one, 
answer "Y" to STRIP PEAK?. If "N", answer 
"Y" to FULL BASE? for option two and, if 
again "N", option three, partial baseline 
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subtraction, is performed. All options 
require very little time and exit to the 
display. 

If, upon entering STRIP the second time, 
the message BAD X POINTS! is issued, the 
user has inverted the left-right X axis 
relationship of the test points. The 
routine immediately cancels the stripping 
operations and exits to the display. 

The user must then reenter STRIP and re¬ 
set the points before proceeding. 

SCALE: 

Scales the CDA to a range of 0-500 
and prints the factors used. 

SMOOTH: 

Smooths the CDA using a slightly modi¬ 
fied SM0T11. SMOT29 is not available. 

INTEGRATE: 

Integrates the CDA after subtracting the 
value of PIVOT LINE from each point. 

The assumed Y axis range is 500, independent 
of the intensity multiplier used to com¬ 
pute LORCAP data. The integral is scaled 
to a range of 500 and the scale down 
multiplier is printed. A modified INTEG 
is used. 

OUTPUT: 

Prints and/or punches the CDA data on the 
device to be specified in answer to TTY 

I/O?. See Section II.G under OUTPUT. 

PLOT: 

Enters the plotter calibration routine. 

See Section II.G under PLOT. 

DIGITAL: 

Enters the digital data input routine to 
allow entry of new data points. See 

Section II.A. 

RESTART: 

Restarts the program. This is equivalent 
to SA=0200. 
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3. Hints on the Use of LORCAP: 

a. A left to right decreasing X axis can be imple¬ 
mented by specifying a negative X RNG and setting 
Xl to the left side value, as usual. This is 
useful in handling NMR data. 

b. An overlapped Y axis may result after using STRIP. 
The normal perspective can be regained by calling 
SCALE. 

c. After startup, the program configuration giving 
peaks or dips and allowing input of equal widths 
can be changed thusly: 

Change C(BSCHNG) to 7041 for dips 

or to 7000 for peaks 

Change C(WIDSWT) to 0 for equal width input 

or to 1 for normal operation 

Restart the program at SA=0201. 

d. When manipulation of the calculated array results 
in its destruction, type 1, call ALTER, and 
enter 0, -1 for peak and parameter numbers to re¬ 
compute the spectrum. This results from the 
fact that the parameters used to compute a spec¬ 
trum are unaffected by all operations except ALTER 
and COMPUTE. 

e. If user was viewing block 0 instead of 1, these 
commands result in the listed symptoms: 

SWAP: Will have no effect whatsoever. 

SUBTRACT: A straight line appears in the 

center of the screen. 

ALTER and COMPUTE: If a spectrum is calcu¬ 
lated, it will go into block 0. 

Although much mention has been made of the rela¬ 
tionship of block 0 and 1, they are, in fact, 
independent of all operations except SWAP and 
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SUBTRACT. Hence, data tapes can be read, and 
computed spectra can be placed into either or 
both arrays. This finds use if it is desired to 
subtract a standard or calibration data tape 
from another data tape, assuming the X axes are 
identical. 

f. The choice for PIVOT LINE for the integrator 
must be fairly accurate, but in the event that 
it was not, the integral can be "corrected" by 
using STRIP to remove a wedge then rescaling, 
if necessary. 

g. To use LORCAP immediately after any one of the 
available MADCAP acquisition routines was exe¬ 
cuted, it is not necessary to have the data 
punched out. It will remain untouched in core 
so long as the LORCAP instruction patch is not 
loaded. 

Load only the first section of LORCAP, use SA=0200, 
answer "N" to CALC'N INPUT? and supply NO. 

POINTS and answer HAVE X DATA? Now, press STOP 
switch on console, LA 0201 and hit START. The 
data will be displayed on the scope and other 
operations can now proceed as usual. 

The author is fully cognizant of the fact that peak 
stripping, extensive smoothing, etc., can be used to 
manufacture data of any desired shape and form; thus a 
note to thesis advisors: Bear a watchful eye on graduate 
students using the program. And, a note to those graduate 
students: Think! Will the results be challenged? Mr. 

John Steichen of Rice University who prompted and aided 
in many of the ideas used herein is acknow]edged for allow¬ 


ing use of these comments. 
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4. Usage of Gaussian Curve Overlay: 

After loading first (and second, if desired) section 
of LORCAP, load the GASCAP overlay. The functional equa¬ 
tion used is changed to 
N 

Yi = Hj * E x P (-(X-Pj)**2/(2* W j**2)) 

j 

with the parameters as described in Section H.l. The over¬ 
lay changes the envelop calculator routine and adds the 
necessary exponentiation routine. Because of the latter, 
the data arrays are cut from 500 points to 464 points; 
i.e., the total data array is 928 points, instead of 1000, 
allowing 464 points each for raw data and calculated data 
or 928 points of either type - all without X data. All 
other operations, starting addresses, etc., remain the 
same as described in Section H.l. 
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III. CONSTRUCTION OF ROUTINES FOR USE WITH MADCAP IV: 

A. General Structure of MADCAP: 

All currently available routines for MADCAP IV 
operate in real time; i.e., program interrupt is not used. 
This allows coding of new routines to be done in a straight¬ 
forward manner without the need to identify and resolve 
concurrent peripheral I/O. Admittedly simple minded and 
inefficient this approach is satisfactory since there is 
little else of a useful nature that could be done simul¬ 
taneously with the few operations that are performed. 

Operating in real time and without recourse to an 
accurately and reproducibly calibrated clock, the averag¬ 
ing routines use memory cycle time to determine real time 
intervals. Hence, for a given machine, reproducibility 
is excellent, but absolute time will vary from one machine 
to another. 

However, if the user desires, analog input routines 
can be written using program interrupt if he constructs 
the necessary interrupt handler and flow controller. It 
must be borne in mind that all I/O thru MADCAP and its 
overlays normally wait for ready or completion flags. 
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The program is centered around the AX08 LAB PERIPHERAL 
which can perform an ADC (analog-to-digital conversion) 
from a prespecified multiplexor channel in about 17 micro¬ 
seconds giving a two's complement, nine-bit signed integer 
result. Nine-bit DAC (digital-to-analog conversion) is 
performed for the scope and plotter driving signals. Hence, 
for both ADC and DAC, the numerical range is 512^o- This 
fact was highly instrumental in the formation of the algo¬ 
rithms for data acquisition, storage, and display. 

Data storage after all types of input should consist 
of single-precision integers in a range of 0-500-^q. 

These are the numbers used for X - Y displays and printout 
and will be characterized within the program by two para¬ 
meters, "PTSWT" in location 66 and "NPTS", the number of 
points, in location 67. 

For Y-only data input (limited to 1000 points), 

PTSWT=1, for digital input wherein the user is supplying 
X data, PTSWT=0 and the maximum number of data points is 
500. The X array is scaled and stored as 0-500 integers 
following the Y array. 

The acquisition and storage of Y data by the analog 
input overlays differ from each other as follows. The 
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CURVE FOLLOWER routine makes an ADC for each point to which 
is added 256^q, making all results positive. That value 
between 0-512 is then stored as single precision. The 
MASS SPEC routine will take 2048 ADC samples adding 256 
to each and store the intermediate sum as a double-precision 
integer. After the sum is complete, the result is divided 
to a 0-512 integer for storage. 

The SIGNAL AVERAGER routine, similarly, takes a 
number of samples (see Table 2) for each data point during 
each scan, adding 256 to each sample, and the samples are 
then summed as a double-precision integer. After sampling 
is complete for that point, the intermediate sum is 
divided to a range of 0-512 and added to the current scan 
sum for that point using a pseudo-18 bit storage scheme; 
i.e., the scan sum uses sesqui-precision storage requir¬ 
ing only 50 per cent more registers than there are data 
points in the scan. 

When the input is terminated, the "18-bit" integers 
are divided by the number of scans taken to give a result 
between 0-512 for storage. 

Note that for analog input which covered the full 
range of the converter, the data must be scaled down from 
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0-512 to 0-500 to prevent wrap-around during scope display 
and plotting. For axes marking, 12 display units are 
reserved on each axis to prevent overdrawing the axes 
onto data. 

B. Some Subroutines Available in MADCAP: 

The more useful and commonly used suboutines are 
addressed indirectly through page zero (locations 20-37). 

A description of each follows: 

(AC=accumulator; FAC=floating point AC in locations 44-46) 
FLOATR (JMS I FLOTER): 

Converts contents of AC to a floating point 
number in FAC. Previous contents of FAC are lost. 

FIXR (JMS I FIXER): 

Converts FAC to an integer leaving result in 
AC and location 45. Values > 2047 or <1-2047 are then 
truncated to those limits. 

READER (JMS I READXY): 

Accepts numeric input (ignoring all but car¬ 
riage returns). It can be implemented in or outside the 
interpreter. In either case, the floating point result is 
in the FAC. If outside, execute JMS I READXY and if in¬ 


side, execute READ. 



STALL (JMS I HOLD): 


Is a variable time delay subroutine using the 
RC clock of AX08. The variable SET, if 0, runs clock 
for 32 pulses, and if non-zero for 8 pulses; i.e., four 
times faster. 

ASRRED (JMS I TELRED): 

Accepts and prints an 8-bit ASCII character 
on the Teletype, leaving the character in the AC. 

INITIZ (JMS I INIT); 

Is called to initialize array pointers and counters 
for data storage. Three anto-index registers are set: 
YIND=YONE (a page zero constant), 

XIND=YONE+NPTS (# of points), 

ZIND=YONE; 

It clears the X display coordinate sum, TXSM; and sets 
a counter, 

CNTR=—NPTS. 

ZIND is used for restorage when a value is extracted 
by YIND and is to be replaced after some modification. 

RETLF (JMS I CRTLFD): 

Issues a carriage-return and line feed. 
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ASKER (JMS I QUERY): 

Allows input of "Y" or "N" in answer to a 
question. All other characters generate "?" and ASKER 
will await proper response. Upon receipt of "Y", it 
exits leaving AC=0. For "N", AC=1. The calling program 
can then branch according to the response. ASRRED is 
used for ASCII I/O. 

UPPEN and DNPEN (JMS I PENUP or PENDN): 

Use enable register and external register of 
AX08 to, respectively, set S0 and clear it to operate a 
relay which controls plotter pen attitude. Each subrou¬ 
tine calls STALL before and after the operation. 

HEDERl and HEDER2 (JMS I HEADRl or HEADR2): 

Output stripped ASCII as generated by MACRO-8 
Assembler. All headings and titles are handled using 
these. There is only one special character, the up-arrow 
(336g), which is used to generate a carriage return-line feed 
in the middle of a text string. Each separate message 
begins with a carriage return-line feed. 

HEDERl takes the content of PC+1 in core as 
the starting address of the message thus allowing the 
message to be stored anywhere in core. If multiple and 
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always consecutive headings are to be printed, the first 
message is printed using JMS I HEADRl followed by address 
and consecutive messages thru JMS I HEADR2 without need 
for message address. For example. 


JMS I HEADRl 
HD I 


JMS I HEADR2 


HDI, TEXT / THIS IS FIRST.t/ 

TEXT/f THEN THIS.'M/ 

When assembled by MACRO-8 and executed will produce 
(CL) 

THIS IS FIRST. (CL) 

(CL) 

(CL) 

THEN THIS. (CL) 

(CL) 
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("CL"=Carriage return-line feed) 

After a call thru HEADRl to initialize the message 
address, text can be printed while in interpreter mode 
by executing HEDIT; e.g., 

• 

JMS I HEADRl 
HDI 


JMS I 7 

HEDIT 

FEXT 

If the messages are not always printed consecutively, 
each must be output by JMS I HEADRl followed by address. 
ADCONV (JMS I ADCON): 

Is a general-purpose subroutine to perform an 
ADC on the channel specified in the AC on entry. The 


result is left in the AC. 



RESET (JMS I DISMOV): 


Resets DAC levels for the scope display to 
contents of XDIS and YDIS which have been set up else¬ 
where. 

PLTINC (JMS I PLTMOV): 

Sets up scope or plotter display depending on 
C(AC). If C(AC)——1 use RESET for scope display; if zero, 
call an incrementer routine to move along the best straight 
line between current position and XDIS, YDIS; and if +1, 
use RESET but keep pen up for a while to allow stabiliza¬ 
tion at new position, then lower the pen. 

ASKIO (JMS I IOASK): 

Requests if user wants TTY I/O. The answer 
sets IOSWT to 0 for Y and to +1 for N. IOSWT is checked 
by both the input and output ASCII routines to choose 
the device before transferring each character. A value 
of 0 uses Teletype and +1 uses high-speed reader/punch. 

Also available are several operations execut¬ 
able while within the floating point package interpreter: 

HEDIT: (See above) 

READ: Input a floating point number ignoring 

non-numeric characters (decimal point 



optional, E format acceptable) 


It is 


logically equivalent to JMS I READXY. 

OUTPUT: Print F. P. number under control of C (62). 

If C(62)=0, use exponential format and 
if non-zero, print integer of that many 
digits. 

NEGATE: Will negate C(FAC). 

SQUARE and SQROOT are available, as usual. 

FLOT: Is used as a memory reference instruction 

to float contents of a single-precision 
operand into the FAC. It supplants FNOR 
as a MRI and FNOR is reduced to the 
status of NEGATE, SQUARE, etc., which do 
not have operands and work only on the 
current C(FAC). 

With the exception of FLOT but including FNOR, these 
interpreter instructions can only be called to a depth of 
one. Hence, a new user function to be called thru the 
interpreter cannot use them but must call those operations 
thru indirect subroutine calls while outside interpreter 


mode. These addresses are as follows: 
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HEDIT: 

JMS 

I HEADR2 

OUTPUT: 

JMS 

I 6 

READ: 

JMS 

I READXY 

NEGATE: 

JMS 

(6000) 

SQUARE: 

JMS 

(6163) 

SQROOT: 

JMS 

(6656) 

FNOR: 

JMS 

(6600) 


With the exception of ASRRED and KBDNTR routines (in 
Basic Package), SERVIS (in LORCAP) and MULT (in Signal 
Averager), all ASCII I/O is handled using routines common 
to the Modified Fl. Pt. Pkg. #2; specifically, output is 
implemented with JMS (7344) and input thru READER using 
JMS (7400). These can be reached thru page zero by JMS I 
OUT and JMS I READXY or READ as described above. Note: 
After a SR start sequence, a TLS is not necessary. 

C. General Layout for New Routines: 

The following logical considerations should be met 
for the writing of new analog input routines. See enclosed 
core map for available space. 

1. Assemble starting address of routine into loca¬ 
tion NOTDIG and into ADRTAB+5, 
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2. Print headings and read in any necessary run 
parameters, 

3. Set NPTS and PTSWT to characterize data; e.g., 
use subroutine DATYP, and any other needed 
variables, 

4. Initialize array storage (JMS I INIT),, 

5. Perform data acceptance, perhaps using JMS I ADCON, 

6. When complete, reduce data such that all Y 
values are = or C500, and 

7. Execute JMP I REETRN to go to display routines. 

For more detailed, yet simple, examples, the 
reader should check the listings for the mass 
spec and line follower routines. 

The basic layout of data treatment routines is similar. 

1. Set up entry logic, one of two ways. 

a. Put starting address of routine in location 
SPOVRl or SP0VR2 and use SR SA=0202 or 0203, 
resp. 

b. or use keyboard command by putting desired 
first letter code character and starting 
address in the tables CHRTAB+X and ADRTAB+X, 
resp., where X is 6, or greater if SMOT 11 or 


29 and INTEG are desired (see below). 
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In either case, the routine is entered by an 
indirect jump. 

2. Initialize arrays, print headings, get input, 
etc. 

3. Perform manipulations. 

4. JMP I REETRN is commonly the last instruction. 

See listing of SMOT29 or INTEG for more detail on 

applications. SM0T11 or 29 use the character "F" in 
location CHRTAB+6 and has its address in ADRTAB+6. INTEG 
uses the character "I" in CHRTAB+7 with the address in 
ADRTAB+7. New routines cannot use the characters Z, S, 

R, P, 0, A, F, or I as the first letter and the address 
must be in a corresponding location in ADRTAB. Z is used 
as a dummy character to fill out the table. Table size 
allows up to five routines to be added to the Basic Package 


in this manner 
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APPENDIX 1: 

INSTRUCTIONS FOR ASSEMBLING MADCAP IV: 

CORE PACKAGE; 

This consists of those locations as shown on the 
accompanying core map. The source tapes do not include 
any overlays or the floating point package. 

The MACRO-8 ASSEMBLER is needed to handle the strip¬ 
ping of the ASCII text. After loading the assembler, the 
NEW ASSEMBLER SYMBOL TABLE is loaded to include AX08 
microinstructions and delete un-needed pneumonics. The 
source tapes are then loaded for PASS I. Only the source 
tapes are re-entered for PASS II (and III, if desired) 

-6 hours are needed to print out PASS III for the entire 
set of MADCAP and overlays onto a Teletype. 

Refer to the MACRO-8 PROGRAMMING MANUAL, DEC-08-CMAA-D, 
for complete instructions in the use of the assembler. 

The high-speed assembler uses the high-speed reader 
for all input and can be modified to use the ASR-33 for 
all output, including the binary tape, by making these 


changes: 
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Change C(603) to 6046 
Change C(603) to 6041 
Change C(3026) to 6044 
Change C(3027) to 6041 
OVERLAYS: 

These are all assembled separately from the CORE 
PACKAGE by means of the USER SYMBOL TABLE, which defines 
the symbols used in the CORE PACKAGE such that overlays 
can be written using the same pnuemonics. 

As above, the MACRO-8 ASSEMBLER is loaded as is the 
NEW ASSEMBLER SYMBOL TABLE. Then the USER SYMBOL TABLE 
is loaded (It should be considered part of the source 
program.), followed by the ASCII source tape(s). The pass 
is terminated with a short tape containing a dollar sign. 

The MADCAP IV BASIC PACKAGE can be made by loading 
the MODIFIED FL. PT. PKG. #2 then the CORE PACKAGE and, 
using ODT-8 (DEC-08-COCO-D-LOW VERSION), punching them 
out with a single checksum. 

The user can also load the various overlays he de¬ 
sires after the CORE PACKAGE and include them in the final 


package. The LORCAP PACKAGE can be made in this manner. 
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APPENDIX 2 t 

MODIFYING MADCAP FOR USE WITH OTHER 

ANALOG I/O DEVICES: 

In any case, data or display variables must not fall 
outside the range of + and -1023 ± 0 . Otherwise, the 
integer arithmetic ingrained in the entire structure 
of MADCAP will give strange results. 

The most likely modification desired, obviously, is 
to live without the AX08 Lab Peripheral. This can be 
done by changing the various routines as described below. 
One common difficulty will be handling ADC’s and DAC's 
with more than 9 bits. Four registers on page zero govern 
display size. These are: XLIM, the left extreme; YLIM, 
the bottom limit; XRNG, which is X axis data range; and 
YRNG, the Y axis range. For the AX08 they were set to 
6, -250, 500, and 500, resp., when assembled. The axes 
marking routine, MARKER, used quadrant fractions of the 
ranges and requires the 6-unit buffer at both ends of each 
axis to draw the enclosing frame. The first page of the 
Core Package listing describes the assembly variables. 

The Core Package can be reassembled with the desired 
values; e.g., double the values to implement 10-bit DAC. 
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If the DAC is 10 bits and the above alterations are 
made, no change should be necessary to use a 10-bit ADC. 

The Signal Averager is limited to 250 scans, however. 

In the display routine, a register, GTX+3, can be 
loaded with CLL RAR, for example, to implement a system 
with a 10-bit ADC but 9-bit DAC. 

Although effort was made to code MADCAP to allow 
these modifications, the author cannot guarantee complete 
success since it was not possible to test them. He is 
willing, however, to confer with interested parties hav¬ 
ing difficulty implementing such changes. 

Here is a breakdown of possibly necessary changes: 

A. Basic Package: 

All analog I/O is handled by five subroutines start¬ 
ing from location 4307. These are RESET, DNPEN, UPPEN, 
STALL, and ADCONV (not used by the Basic Package it¬ 
self) and are described in Section III. It is fairly 
simple to overlay the necessary IOT's for all except 
perhaps those implementing the enable and external regis¬ 
ters and the RC Clock of the AX08. STALL can be 
emulated with a nested pair of ISZ loops or if an 
analog input is available; i.e., a potentiometer 
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delivering + to -1 volt or so into an ADC, it can be 
used to get variable delay times. The delay period 
should be in the range of 0.1 to 2 seconds. 

Although quite awkward, UPPEN and DNPEN can be 
written to wait for a Teletype flag to allow manual 
lifting and lowering of the pen. 

B. LORCAP: 


C. 



D. 


Use is made of all the analog I/O routines and 
four analog inputs going into consecutive multiplexor 
channels. The first channel is written in as 34g but 
this can be reset by changing C(PL34) to the desired 
channel. 

Curve Follower Input: 

It uses only those routines in the Basic Package. 
Mass Spec Input: 

The external register is used to sense a Schmitt 
trigger. This can be implemented by monitoring an 
ADC channel for any desired preset level as detected 
by the software. Change locations L00P3, L00P3+1, 


and PTLOOP + 7 to 11. 


E. Signal Averager: 

All special analog I/O is performed in sub- 
routine AVGMR. What is needed here is a starting 
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trigger mechanism, another programmed but accurate 
delay loop, and, if desired, a means for the com¬ 
puter to signal the end of a scan. The changes here 
are reasonably easy to make with a variety of hard¬ 
ware approaches. 
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APPENDIX 3 s 


Functional Flow Chart 

for Basic Package; 
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APPENDIX 4: 


CORE MAP FOR MADCAP ROUTINES: 

(All locations are in octal.) 


Basic Package: 


2-144 : 
200-377 : 
2400-2577: 
4000-5377: 
5400-7577: 


FPP registers, subroutine references, etc. 
Main calling routine and display 
MODIFIER routine 
Operational routines 
Modified FPP 


The following assume all registers of the Basic Package: 
LORCAP Package: 


145-177 : Subroutine reference, constants, etc. 
400-2350: Reserved for data 
2350-2462: Reserved for peak parameters, contains 
initial dialog routine on startup 
2463-3777: Operational routines 

Note that the MODIFIER is lost. 


Averager: 

145-177 : Subroutine references, constants, etc. 
3400-3777: Operational routines 

Also, the code used by the Basic Package for digital input 
and Y axis tic marks is overlayed. MODIFIER is destroyed 
by a 1000-point run. 
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Mass Spec: 

3600-3777: Operational routines 

4160-4167: Stripped ASCII text 

Curve Follower: 

3600-3777: Operational routines 

SMOT11 and 29: 

3200-3322 and 3355, resp. 

INTEG: 

3000-3144 
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APPENDIX 5: 

Samples of Program Output: 

(Text in parentheses is added comments; underlined text is 
user input) 

(Basic Package, SM0T11, and INTEG were loaded with BIN 
Loader. Then used SA=0200.) 


MADCAP IVI 

TITLE: EXAMPLE OF USE OF 

BASIC PACKAGE 

ANALOG INPUT7N 

NO. POINTS= 100 
HAVE X DATA? IT 
PTS. TO SKIP = 0 

RANGE & MIN 
Y: + 500 - 0 


DATA: 

TTY l/0?N 


S: 

MIN= 46 
MAX= 335 
MPY= 1.730 
P: 

CALIBRATE: 
LINE PL0T7Y 

YXO 7 

_7_ ? _?OM 

G 


(Terminated with ALT 
MODE) 

(Want data tape input) 


(Spaces between + or - 
sign and number are 
ignored) 


(Use high speed reader) 
(Scope display began after 
input) 

(Scale the data) 


(Plot the data) 


(Connect points with best 
straight line) 

('?' from illegal char¬ 
acters - spaces) 

('M' to draw axes and 
tic marks) 

('G* begins plot) 
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APPENDIX 5: continued 


TURN PLTR OFF! 


0 ; 

TTY l/0?N 


R: 


NO. POINTS « 100 
HAVE X DATA7N 
PTS. TO SKIP = 0 


DATA: 

TTY l/0?N 

F: 

I: 

PIVOT LINE (0-50jZhl0 
MPY=0.233296E-01 


(Plot is complete? scope 
display restarts) 

(Output the data) 

(Use high speed punch) 

(Re-enter data from tape) 


(Use high speed reader) 

(Filter the data) 
(Integrate it) 


(Scale down multiplier) 


P: 


(Plot data) 


CALIBRATE: 
LINE PL0T7N 
G 


(Make point plot) 


TURN PLTR OFF! 
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APPENDIX 5; continued 


(Loaded LORCAP IV, first section, then used SA=0200) 
LORCAP IV! 


TITLE: SAMPLE OF LORCAP IV 

INITIAL DATA BLOCK NO. = 0 
WANT ABSORPTION 'DIPS 1 ?N 
ALL WIDTHS EQUAL7Y 

CALC'N INPUT7N 

NO. POINTS^ 100 
HAVE X DATA7N 
PTS. TO SKIP - 0 

RANGE & MIN 
Y: 5gTgT 0 

DATA: 

TTY I/07N 

SCALE: 

MIN= 46 
MAX= 335 
MPY= 1.730 
OUTPUT: 


TTY I/07Y 




2 

2 

3 

5 

5 

14 

14 

16 

17 

21 

36 

40 

45 

48 

54 

95 

106 

116 

130 

144 

279 

310 

343 

375 

410 

490 

472 

446 

415 

382 

202 

182 

163 

147 

131 

73 

66 

61 

55 

50 

28 

26 

22 

21 

19 

9 

9 

7 

5 

5 


(Enter data in BLOCK 0 ) 
(Use 'peak' spectra) 


(Input the data) 


(Use high speed reader) 
(Scale the block 0 data) 


(Output block 0 data) 
(Use Teletype) 


7 

7 

9 

10 

12 

22 

24 

28 

29 

33 

59 

64 

71 

78 

87 

159 

178 

197 

220 

246 

441 

467 

488 

498 

500 

349 

317 

285 

256 

230 

119 

107 

97 

88 

80 

45 

42 

38 

35 

31 

17 

16 

14 

12 

10 

3 

3 

2 

2 

0 










APPENDIX 5: continued 


STRIP: 

(Put test points on scope) 

ST: 

(After positioning them. 


re-enter STRIP) 

STRIP PEAK7N 

FULL BASE7Y 

(Restored sloping baseline) 

SCP? 

(Hit a RUBOUT) 

SCALE: 

(Scale Block 0 data) 

MIN= 229 

MAX= 380 

MPY= 0.821 

ST: 


STRIP: 

(Enter STRIP again) 

STRIP PEAK7N 

FULL BASE7N 

(Remove portion of baseline) 

SCALE: 

(Scale data) 

MIN“ 0 

MAX 55 500 

MPY“ 1.000 

FILTER:? 

(Illegal command) 

SMOOTH: 

(Smooth data) 

SM: 

(Smooth data again) 

ST: 

(Enter STRIP) 

ST: 

STRIP PEAK7N 

FULL BASE7N 

(Remove more baseline) 

SMOO: 

(Smooth again) 

(NOTE: 

All of the above was performed 


on BLOCK 0.) 


01CO : (Typed 1 and entered COMPUTE) 


NEW CALCULATION: 

NO. POINTS =100 (Must be same as Block 0) 

HAVE X DATA7N 
NO. PEAKS=2 ' 
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APPENDIX 5: continued 


X RNG, XI, INT.MPLR, BASE: 
10 10 5 10 

HGT,H-WIDTH,POS: 

80 1 15 

70 17 


COOL IT. 
0101011SC: 


MIN= 65 

MAX= 449 

MPY= 1.302 
01011SUBTRAAAAAACT : 

SU: 


ALTER: 

PEAK & PARA #’S : 0 2 
TO: 15 

PEAK & PARA # 1 S: 1 1 
TO: 85 

PEAK & PARA # 1 S: 0 -1 

COOL IT. 

AL: 

PEAK & PARA # 1 S: 0 0 


(No width needed, ’1 * is 
assumed) 


(Alternately displayed arrays 
then SCALED Block 1) 


(Typed 1 then, effectively, 
SU: for SUBTRACT:) 
(SUBTRACT again to restore 
array) 

(ALTER para's) 

(Alter base) 

(Alter height of peak 1) 

(Exit ALTER) 

(Computing Spectrum) 

(Enter ALTER again) 

(Get para listing) 


X RNG, XI, INT.MPLR, BASE: 

0.100000E+02 0.100000E+0? 0.500000E+01 0.150000E+O? 

HGT,H-WIDTH, POS: 

0.850000E+02 0.100000E=02 0.150000E+02 
0.700000E+02 0. 100000E+01 0.170000E+02 

AL: (Enter again) 

PEAK & PARA #'S: 1 0 (List only first peak Paras) 

0.850000E+02 0.100000E+01 0.1500000E+02 


AL: 

PEAK & PARA # 1 S: 1 2 
TO: 1.5 

PEAK & PARA #'S: 1 3 
TO:16 

PEAK & PARA #'S:2 1 


(Enter again) 

(Width of peak 1) 

(Position of peak 1) 


(Height of peak 2) 
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APPENDIX 5: continued 


TO: 75 

PEAK & PARA # 1 S: 1 4 
PEAK & PARA # 1 S: - 1 
PEAK & PARA #'S:25 -1 


COOL IT. 

ST: 

ST: 

STRIP PEAK7N 
FULL BASE7Y 

jglSM? 

SC: 


MIN= 13 
MAX= 468 
MPY= 1.099 
IN: 

PIVOT LINE (0-500)=0 
MPY~ 0. 125780E-01 

101AL: 

PEAK & PARA #'S: 0 -1 

COOL IT: 

SU: 

SC: 


MIN= - 42 
MAX= 343 
MPY= 1.299 
IN: 


PIVOT LINE (0-500 ) =250 
MPY= 0.490629E-01 
RESTART: 

LORCAP IV! 


TITLE: 

CALC'N INPUT7N 

NO. POINTS^ 100 
HAVE X DATA7N 
PTS. TO SKIP - 10 


(PARA 4 is illegal) 
(PEAK-1 is illegal) 

(PEAK 25 is ignored; PARA 
-1 for exit) 


(Enter STRIP) 


(Remove sloping baseline) 
(Smooth Block 1) 

(Scale it) 


(Integrate Block 1) 


(Enter ALTER) 
(Recompute spectrum) 


(Get residual 
(Scale them) 


s) 



(Integrate residuals) 


(Restart program) 


(Hit ALTMODE) 


(Skip first 10 points on tape) 
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APPENDIX 5: continued 


1SWAP 

P; ? 

0PL; 


(Invalid command) 
(Plot Block 0) 


CALIBRATE: 

LINE PL0T7Y 

QXQXQXQXOXOYOYOYOYOYQYOYOM 


(Calibrate plotter, draw 
axes and tics) 


G 


TURN PLTR OFF! 


01PL: 


(Plot Block 1) 


CALIBRATE: 

LINE PLOT7N 
G 

TURN PLTR OFF! 

WELL, THAT IS ALL THERE IS TO IT!!!:? 

(This was typed during scope display and '?' means 
it was an invalid command, obviously.) 
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APPENDIX 5: continued 


(Loaded first section of LORCAP, then Gaussion overlay. 
Used SA=0200) 


GASCAP IV! 

TITLE: SAMPLE OF USE OF GAUSSIAN OVERLAY 

CALC 1 N INPUT7N 

NO. POINTS = 100 
HAVE X DATA7N 
PTS. TO SKIP = 10!_0 

RANGE & MIN 

Y: 500 0 

DATA: 

TTY l/o?N 

y 

CO: 

NEW CALCULATION: 


(Wanted 0, not 10; hit RUBOUT) 


NO. POINTS = 100 
HAVE X DATA7N 
NO. PEAKS = 1 

X RNG, XI, INT.MPLR, BASE: 
100 0 5 0 

HGT,H-WIDTH, POS: 

90 6 50! _55 


(Wanted 55, 
RUBOUT) 


not 50; hit 


COOL IT. 
011000DI: 


(Computing spectrum) 

(Get new data into Block 0) 


NO. POINTS * 100 
HAVE X DATA7N 
PTS. TO SKIP = 0 
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APPENDIX 5 


RANGE & MIN 

Y:500_0 

DATA: 


continued 


TTY l/0?N 


(All operations same as 
LORCAP) 
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APPENDIX 5; continued 

AVERAGER! 

TITLE: SAMPLE OF SIGNAL AVERAGER 

(Loaded Basic Package and Averager overlay; used 
SA-0200) 


MPX CHAN-28 

PT. INDEX = _1 

SEC/SCAN=.6 

PT. INDEX-1. 

SEC/SCAN- 1.0 

MAX SCANS = 3 

DELAY (SEC)= 0 

(Must be decimal) 

(50 points) 

(Too little time) 

• 

(Input commenced and three 
scans were taken) 

A: 

PT. INDEX =1 

SEC/SCAN=10 

MAX SCANS =2 

DELAY (SEC) =0_ 

(Restart analog input) 

(Input began on first trigger) 

NO. SCANS= 1 

(Bit 0 of SR was set to 
interrupt) 

dddmmmmmmddds 

('D' to decrease size,'M" 
to magnify, 'S' to stop) 

MORE SCANS ?Y 

(Continued averaging, went 
to completion and scope disF^y 
began.) 

I: 

F: 

PI: 

RE: 

PT. INDEX=3 

SEC/SCAN=60 

MAX SCANS=_1 

DELAY(SEC)=0 

(No action; INTEG not loaded) 

(No action; FILTER not loaded) 
(Ditto: invalid command) 

(Restart, input routine) 

(500 points) 


Remaining operations as in Basic Package 
(Reloaded Averager overlay only; used 
SA=0400 to reset machine cycle time.) 

MACHINE CYCLE TIME (SEC) = 1.5E-6 


(Program continues as if SA=0200) 










4 . 


OVERLAY PLOT OF EIGHT PEAK 
LORCAP ON DATA POINTS,WITH 
RESIDUALS (SMOOTHED): 









DECtape ^1 contains copies of the source and binary tapes. 

M4BIN. 1 is equivalent to tape 1 of 4 
M4BIN.2 is equivalent to tape 2 of 4 

and 

M4SRC. 1 is equivalent to tape 1 of 6 
M4SRC.2 is equivalent to tape 2 of 6 

etc. 

To get punch copies from the DECtape using PIP on the PDP-10, the source files 
can be directly transferred, e.g., 

PTP: 4— M4SRC.1 

but the binary files must always be tranferred using image mode, e.g., 

PTP:4— M4BIN. l/l 
and also, 

D$K;/X4—M4BIN. 1/1 
etc. 

DECtape ^2 is an edited, complete PAL-10 listing of all the program segments. A 
printer listing is acquired thusly, using PIP: 

LPT: 4— DTAN:M4LST 





